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Abstract 

Vortices in non-Abelian gauge field theory play important roles in con- 
finement mechanism and are governed by systems of nonlinear elliptic 
equations of complicated structures. In this paper, we present a series 
of existence and uniqueness theorems for multiple vortex solutions of the 
non-Abelian BPS vortex equations over M 2 and on a doubly periodic do- 
main. Our methods are based on calculus of variations and a fixed-point 
argument. The necessary and sufficient conditions for the existence of a 
unique solution in the doubly periodic situation are explicitly expressed in 
terms of several physical parameters involved. 



1 Introduction 



Vortices have important applications in many fundamental areas of physics. 
For example, in particle physics, vortices allow one to generate dually (electri- 
cally and magnetically) charged vortex-like solitons [17] [29] [32] known as dyons 
[31] [34] [35]; in cosmology, vortices generate topological defects know as cosmic 
strings [13] [28] which give rise to useful mechanisms for matter formation in the 
early universe. Besides, both electrically and magnetically charged vortices arise 
in a wide range of areas in condensed-matter physics including high-temperature 
superconductivity [9] [26], optics [6] [27] [30], and so on. 

Mathematically, Chern-Simons theories in (2+l)-dimensions are introduced 
to accommodate electricity. The equations of motions of various Chern-Simons 
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vortex moldels are hard to approach even in the radially symmetric static cases. 
However, since the discovery of the self-dual structure in the Abelian Chern- 
Simons vortex model [14] [15] in 1990, there came a burst of fruitful works 
on Chern-Simons vortex equations, non-relativistic and relativistic, Abelian and 
non- Abelian [7] [8]. For example, Aldrovandi and Schaposnik [3] [24] found the 
non-Abelian vortex solutions when gauge field dynamics is solely governed by a 
Chern-Simons action and the symmetry breaking potential is six-order in order to 
ensure self-duality and supersymmetric extension, in the presence of a set of ori- 
entational collective coordinates. Furthermore, the existence of Chern-Simons- 
Higgs vortex solutions was proved in (2+l)-dimensions with internal collective co- 
ordinates [18]. The existence of topological solutions for relativistic non-Abelian 
Chern-Simons equations involving two Higgs particles and two gauge fields was 
proved through studying the full M 2 limit of a coupled system of two nonlin- 
ear elliptic equations [21]. In 2008, Aharony, Bergman, Jafferis and Maldacena 
developed the so-called ABJM theory [2] in terms of three dimensional Chern- 
Simons-matter theories with gauge groups U(N) x U(N) and SU(N) x SU(N) 
which have explicit J\f—6 superconformal symmetry. Before long, Auzzi and Ku- 
mar [5] find half-BPS vortex solitons, at both weak and strong couplings, in this 
theory. 

More recently, the existence of solutions for Abelian Chern-Simons equations 
involving two Higgs particles and two gauge fields on a torus was proved by Lin 
and Prajapat [20]. Using the methods of monotone iterations, a priori estimates, 
degree-theory argument and constrained minimization, multiple vortex equations 
in U(N) and SO(2N) theories were discussed [12] [22] [23] and a series of sharp 
existence and uniqueness theorems were established. Lieb and Yang [19] discussed 
non-Abelian vortices in supersymmetric gauge field theory, over doubly periodic 
domains, via a highly efficient direct minimization approach. These studies unveil 
a broad spectrum of systems of elliptic equations with exponential nonlinearities 
and rich properties and structures, which present new challenges. 

In this paper we will concentrate on the non-Abelian BPS vortex equations 
derived by Auzzi and Kumar [5] in a supersymmetric Chern-Simons-Higgs the- 
ory formulated by Aharony, Bergman, Jafferis and Maldacena [2], known as the 
ABJM model. In terms of the methods of [19] [22] [23] [33], we obtain the exis- 
tence and uniqueness of a multiple vortex solution. Meanwhile, it is hopeful that 
our method may be explored further to study various multiple vortex equations, 
arising in non-Abelian gauge field theory, of more diffcult structures. 

The content of the rest of paper is outlined as follows. In Section 3, we prove 
the existence and uniqueness of a multiple vortex solution realizing an arbitrarily 
prescribed vortex distribution over M 2 , applying the variational method of Jaffe 
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and Taubes [16] used for the Abelian Higgs model. In Section 4, we prove the 
existence of a multiple vortex solution over a doubly periodic domain under a 
necessary and sufficient condition explicitly stated in terms of some physical cou- 
pling parameters, by a multi-constrained variational approach. In Section 5, we 
prove the existence of multiple vortex solutions in a doubly periodic domain by 
a fixed-point method, where we apply the technique of the maximum principle 
and the Poincare inequality. Furthermore, in Section 6, our methods are shown 
to be equally effective in treating the existence and uniqueness problems for the 
multiple vortex solution induced from independently prescribed distributions of 
zeros of two complex scalar fields, instead of one. 



2 Non- Abelian vortex equations 



Recently developed by Aharony, Bergman, Jafferis and Maldacena, known 
also as the AB JM model [2] , is a Chern-Simons theory within which the matter 
fields are four complex scalars, 

C 1 = (Q 1 ,Q 2 ,R 1 ,R 2 ), 7 = 1,2,3,4, (2.1) 

in the bifundamental matter field (N, N) representation of the gauge group 
U(N) x U(N), which hosts two gauge fields, and B^. The Chern-Simons 
action associated to the two gauge group and of levels +k and — k is given 
by the Lagrangian density 



Cos = — e^Tr (A^A, + -A,A U A, - B lx d v B 1 - -B,B U B^ , (2.2) 

where the gauge-covariant derivatives on the bifundamental fields are defined as 

D^C 1 = d^C 1 + lA^C 1 - iC l B^ 7=1,2,3,4. (2.3) 

The scalar potential of the mass deformed theory can be written in a compact 
way as [11] 

V = Tr(M at M a + N at N a ), (2.4) 



where 



M a = pQ a + —(2Q [a QlQ /3] + RpR ] pQ a - Q a B} p R p 

+2Q?R\R a - 2R a R\Q?) , (2.5) 

N a = - P R a + ^-(2R [a RlR p] + Q p Q\R a - R a Q\Q p 

+2RPQlQ a -2Q a QW), (2.6) 
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where the Kronecker symbol e a/3 (a, /3 — 1, 2) is used to lower or raise indices, and 
p > a massive parameter. Thus, when the spacetime metric is of the signature 
(H ), the total (bosonic) Lagrangian density of ABJM model can be written 

as 

C = -C cs + MID.C^ID^C 1 ]) - V, (2.7) 

which is of a pure Chern-Simons type for the gauge field sector. As in [5], we 
focus on a reduced situation where (say) R a = 0. Then, by virture of (2.5) and 
(2.6), the scalar potential density (2.4) takes the form 

A-TC 

V = Tr(M a ^M a ), M a = P Q a + —(Q a QlQ? -Q?QlQ a ). (2.8) 

The equations of motion of the Lagrangian (2.7) are rather complicated. However, 
in the static limit, Auzzi and Kumar [5] showed that these equations may be 
reduced into the following first-order BPS system of equations 

DqQ 1 - iW 1 = 0, D X Q 2 - iD 2 Q 2 = 0, (2.9) 
D 1 Q l = 0, D 2 Q l =0, D Q 2 = 0, W 2 = 0, (2.10) 

coupled with the Gauss law constraints which are the temporal components of 
the Chern-Simons equations 

A e A-7 F 0A) = i(Q«[£)MQa]t _ [D"Q"]Q°*), (2.11) 
47T 

Air 



where 



F#> = d,A u -d u A^ + i[A^,A u ], 
= d.Bv-d^ + iiB^B,,}, 

W 1 = pQ' + ^iQ'Q^Q'-Q'Q^Q 1 ), 
k 

W 2 = pQ 2 + 2 -^{Q 2 Q l ^Q l -Q l Q^Q 2 ), 



provided that [5] one takes that Q l assumes its vacuum expectation value 

Ql = V S diag (°' " ' v/F ~^' ^ F3T ) > ( 2 - 13 ) 



the non-trivial entries of Q 2 are given by (N — 1) complex scalar fields k and 
(fie (£ — 1, • • • , N — 2) according to 

Q 2 n,n-i = Q 2 N-e,N-e-i = (2-14) 
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and the spatial components of the gauge fields Aj and Bj (j — 1, 2) are expressed 
in terms of (N — 1) real-valued vector potentials a = (aj) and b = (bj) (j = 
1,2; £ = 1, - ■ ■ , N — 2) satisfying 

A j = B j = diag(0,aJ r - 2 ,...,oJ,6 j ), J = 1,2. (2.15) 

We now consider the solution for the N = 3 case. The ansatz for the bifun- 
damental scalars approaching the Higgs vacuum at infinity is 



Q 1 



Q 2 = \ V2k \ , (2.16) 




Aj = B< ; = I a, I , ?' = 1,2. 



Where k is a real-valued scalar field, a complex-valued scalar field, and aj and 
ft, are two real-valued gauge potential vector fields. 

Define cijk = djdk — dk&j and set A = 4p 2 . Then the vortex equations without 
assuming radial symmetry are 

(d 1 +id 2 )K = i(a 1 + ia 2 )n, (2-17) 
(d 1 + id 2 )<f) = -i([ ai + ia 2 ] - [h + ib 2 ])(f). (2.18) 

a 12 = -^-H 2 -!), (2.19) 
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-A(|0| 2 -1). (2.20) 



We shall look for solutions of these equations so that k never vanishes but (p 
vanishes exactly at the finite set of points 

Z={ Pl ,p 2 ,---,p n }. (2.21) 

A solution is called an n- vortex solution [16]. 

To facilitate our computation, it will be convenient to adopt the complexified 
derivatives 

0=1^-^), d= l -{d l + id 2 ), (2.22) 

and the notation 

a = cii + ia 2 , b = b 1 +ib 2 . (2.23) 
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As a consequence, away from Z, the equations (2.17) and (2.18) become 

d\nK=--a, d\n<p = --(a- b), (2.24) 

which allow us to solve for a, b to get 

a = 20 In k, a- 6 = 2i9 In 0. (2.25) 

Using 

a 12 = -i(da-da), (2.26) 
(2.19), (2.20), (2.25), and the fact that dd = dd = ±A, we have 

ai2 = -Aln/t. (2.27) 

Likewise, we have, away from Z, the relation 

612 = ai2 - In |0| 2 = -^A(ln k 2 + In |0| 2 ). (2.28) 

Set u = IriK 2 and v — In |0| 2 and note that |0| behaves like \x — p s \ for x near 
p s (s — 1, • • • , n). We see that u and v satisfy the equations 

Au = X(2e u -e v -l), (2.29) 

n 

Au + Av = 2X(e v - 1) +4tt^5 Ps (x), (2.30) 

s=l 

where we have included our consideration of the zero set Z of (f> as given in 
(2.21). 

3 Solution on full plane 

In this section, we prove the existence and uniqueness of the solution of the 
system of equations (2.29) and (2.30) over M 2 satisfying the boundary condition 

u — > 0, v — > as |x| — > oo. (3.1) 

To proceed further, we introduce the background function [16] 

n 

v (x) = -^ln(l + r\x-p s \~ 2 ), r>0. (3.2) 

s=l 
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Then, we have 

n n 

Av = -h(x)+An^26 p .(x) 1 h(x) =4^ . 

8=1 3 = 1 ^ + |X Psl I 



(3.3) 



Using the substitution v = v + w, we have 



Au = X(2e u -e vo+w - 1), (3.4) 

A(u + w) = 2X(e V0+w -l) + h(x). (3.5) 

Taking / = u + w, we change (3.4) and (3.5) into 

Au = X(2e u - e vo+f ~ u - 1), (3.6) 

A/ = 2X(e V0+f - u -l) + h{x). (3.7) 

It is clear that (3.6) and (3.7) are the Eulur-Lagrange equations of the action 
functional 

I(u,f) = / 2 {^|V«| 2 + ^|V/| 2 + (2(e«-l)- W ) 

+ ( e vo+f-u _ e v ^ + _ -q j} dx 
2X ) 

It is clear that the functional / is a C^-functional for u, f E W 1 ' 2 (M 2 ) and its 
Frechet derivative satisfies 

DI(uJ)(u,f) = l 2 {\\Vu\ 2 + ^\Vf\ 2 + e^(e^-l)(f-u) 

+2 (e« - l)u + (e v ° -!)(/-«) + ^f}dx. (3.9) 



h 

2X 

Since 



|Vw| 2 + |V/| 2 = 2|Vw| 2 + |Vuf + 2(Vu, Vw), (3.10) 

Hence 

|V«| 2 + |V/| 2 < 3|V«| 2 + 2|V«;| 2 < 3(|V«| 2 + |V^| 2 ). (3.11) 
On the other hand, we have 

|Vn| 2 +|V/| 2 > 2\Vu\ 2 +\Vw\ 2 -2\{Vu,Vw)\ 

> {2--)\Vu\ 2 +{l-e)\Vw\ 2 , (3.12) 

e 

for any e E (|, 1). 
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Taking e — |, we get 



\Vu\ 2 + |V/| 2 > i|V«| 2 + ±\Vw\ 2 > |(|W| 2 + |VH 2 ). (3.13) 



Similarly, we have 



\{u 2 + w 2 ) <u 2 + f <3(u 2 + w 2 ). (3.14) 



As a consequence of (3.9), (3.13) and (3.14), we obtain 
£>/(«, /)(«, f)~^J {l V«| 2 + |Vuf 

> / {e^(e w -l)w + 2(e u -l)M+(e^-l)w + ^(w + u0i 

= / {(e^(e w -l)+e ,;o -l + ^)w + (2(e M -l) + 

Jr2 2A 2A J 



iw(e W0+,u - 1 + + u(2(e u - 1) + ^r)\dx 
2a 2a J 

= M 1 (w) + M 2 (u). (3.15) 

As in [16], we decompose w and u into their positive and negative parts, w = 
w + — w_ and u = u + — u_, where q + = max{g,0} and g_ = — min{g, 0} for 
q ei Using the elementary inequality 

e'-l>t, teR, (3.16) 

we have 

e v ° +w -l + ^>v + w + ^, (3.17) 

Mi(w + ) > [ w 2 + dx+ [ w + (v + ^-)dx 
Jr 2 Jr 2 2A 

1 f 1 f h 

> / ^ / (vq+ )2dx _ (318) 

z JR 2 ^ JR 2 ZA 



which leads to 



On the other hand, using the inequality 



l- e -*>-l_, t>0, (3.19) 



we have 



w_(l - ^- - e Vo ~ w -) = w4l-^- + e V0 (l-e~ w -)-e Va ) 
2a 2a 

W - -(l-A) + _?b_(l- e -- *).(3.20) 



l + 2\' 1 + w- 2A' 
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In view of (3.3), we see that we may choose r > large enough so that 

h.(r\ 

< 1, xe R . (3.21) 



A 

Since 1 — e v ° and h both lie in L 2 (R 2 ), we have 



/ —l-e v "- — dx<e — dx + C(e), (3.22) 

Jr 2 1 + w- 2A " U l + w- K h K J 

where e > may be chosen to be arbitrarily small. Combining (3.20)-(3.22), we 
obtain 

If w 2 

M 1 (-w.) > - \ ——dx - d(e), (3.23) 
4 Jr 2 1 + "Un- 
provided that e<\. From (3.18) and (3.23), we get 

If w 2 

M 1 (w)>- —dx-C, (3.24) 

where and in the sequel we use C to denote an irrelevant positive constant. Similar 
estimates may be made for M 2 {u). Thus, (3.15) gives us 

DI(u, /)(«,/) - / {|Vm| 2 + |VH 2 }<fc 
>7/f^ + T^V-C. (3.25) 



4 JR2 VI + |u| 1 + |w. 
We now recall the well-known Gagliardo-nirenberg inequality 

/ u 4 dx <2 [ u 2 dx [ \Vu\ 2 dx, u e W h2 (R 2 ). (3.26) 
Jr 2 Jr 2 Jr 2 

Consequently, we have 

'4' = (/LttL< i+ '"'>i«H 2 



(1 + lul) 2 ^ j R 



- — rr^dx / (1 + \u\) \u\ dx 



r 2 

< 2 



< 4 



< 



/ ~, — —, — ttt: dx / (u 2 + u 4 )dx 
J R2 (l + \u\) 2 J R2 K } 

/ 7 i — TT^dx / u 2 dx( 1 + / |Vu| 2 rfx J 

Jr 2 (i + mr Jr 2 v j R 2 / 

KI^) 2 + c ( 1 + iI(i^^ 4 

+[/ IVwI 2 ^] 4 ). (3.27) 

JR 2 / 



As a result of (3.27), we have 

(I u 2 dx) 2 <c(l+ [ \Vu\ 2 dx + [ - — ^—-r-dx). (3.28) 

Applying (3.28) in (3.25), we arrive at 

DI(u, f)(u, f) > Ci(||w|| w i,2 (M2) + ||w|| w i,2 (M 2)) - C 2 . (3.29) 

Thus, using (3.11), (3.13) and (3.14) in (3.29), we conclude with the coercive 
lower bound 

DI(u, f)(u, f) > CiflMIwi,^) + 11/11^1,2^)) - C 2 . (3.30) 

With (3.30), we can now show that the existence of a critical point of the 
action functional (3.8) follows by using a standard argument as in [33]. 
In fact, from (3.30), we can choose R > large enough such that 

mi{DI(u,f)\u,fe W 1 ' 2 ^ 2 ), |MUi, 2(K2) + WfWw^m = R}>1 (3.31) 

(say). Now consider the minimization problem 

7] = mm{I(u,f)\ ||m||wi.3(r2) + 11/11^1,2^2) < R}. (3.32) 

Let {(uk, fk)} be a minimization sequence of (3.32). Without loss of generality, 
we may assume that {(«&, /&)} weakly converges to an element (u, /) in W 1,2 (M. 2 ). 
The weakly lower semi-continuity of / implies that (u, /) solves (3.32). To show 
that (u, /) is a critical pint of /, it suffices to see that it is an interior point. That 
is, 

IMIwn. 2 (R 2 ) + ll/HwL^R 2 ) < R- (3.33) 

Suppose otherwise that ||u|| w i,2( R2 ) + ||/||iv 1 , 2 (r 2 ) = R- Then for t G (0,1) the 
point (1 — t){u, /) is interior which gives us 

I((l-t)u,(l-t)f)>r] = /(«,/). (3.34) 

On the other hand, we have 

& /((i -mu, /»-/(»,/) = | /((1 _ t)(Bi/))U 

= -DI(uJ)(u,f) <-l. (3.35) 
Consequently, if t > is sufficiently small, (3.34) leads to 

I((l-t)(u, f))<I(u,f) = ?7, (3.36) 
which contradicts (3.34). Therefore, the existence of a critical point of / follows. 
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Note that the part in the integrand of / which does not involve the derivatives 
of u and / may be rewritten as 

Q(u, f) = 2{e u -l)-u + e v ° + f~ u - e"° + ( A _ 1)/, (3.37) 

whose Hessian is easily checked to be positive definite. Thus, the functional / is 
strictly convex. As a consequence, / can have at most one critical point (u, f) in 
the space W li2 (R 2 ). 

To proceed further, we now show that the following claim holds. 

Claim: If g e W 1>2 (R 2 ), then e 9 — 1 E L 2 (R 2 ). 

We first recall the Sobolev embedding inequality in two dimensions [10]: 

\\g\\LkQ&2) < ( 7r (^ — )) 2k WaWw^m, k>2. (3.38) 
On the other hand, the MacLaurin series leads to 



00 ok 



(e*-l)W + £A^A (3.39) 

k=3 

Combining the above with (3.38), we have, formally, 

2 2 ^2 k -2/ k-2\^ k 

\\ e - i|Il2(r 2 ) < ||5lLfc(R2) + — — y — 2 — ) WdWw 1 ' 2 ^ 2 )- (3.40) 

fc=3 

Setting 

2 k -2f k-2\^ k 



and applying the Stirling formula, 

k\ ~ V2^k k+ h~ k (k -> oo), (3.41) 

we have 



^W^2 ( k-2\Y 2 

- k 'a k ~ 



y z^ — z / K — z\-*k 

r U — — \\9 Ww 1 

ke- 1 (2kn)^y 2 J 

/ k — 2\ h 

~ 2e^\\g\\ w ^ m {—) ^0<1 (fc-xjo). (3.42) 

Thus we have shown that (3.40) is a convergent series, which verifies our claim. 

We now continue our work. Noting vo,h G L 2 (R 2 ) and using the claim, we 
see that the right-hand side of (3.6) and (3.7) belong to L 2 (IR 2 ). We may now 
apply the standard elliptic theory to (3.4) and (3.5) to infer that u,w e W 2 ' 2 (M 2 ). 
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In particular, u, w and |Vw|, \ Vw\ approach zero as \x\ — > oo, which renders the 
validity of the boundary condition (3.1). 

Finally, we derive the decay rates for u,v and |Vw|, |Vi>|. Consider (2.29) and 

(2.30) outside the disk D R = jz G R 2 \x\ < where 

R > max| \p s \ s — 1, 2, • • • , n j. 

We rewrite (2.29) and (2.30) in R 2 \ D R as 

Au = X{2e u -e v -l), (3.43) 
Av = X(-2e u + 3e v - 1). (3.44) 

By computation, we have 

A(u 2 + v 2 ) = 2{\Vu\ 2 +\Vv\ 2 )+4Xu{e u -l) + 6Xv{e v - 1) 

-2Xu(e v - 1) -AXv(e u - 1), xGl 2 \D«. (3.45) 

Noting -u, v — )> as |x| — > oo, for any e : < e < 1, we can find a suitably large 
R e > R so that 

A(w 2 + w 2 ) > 4Aw 2 + 6Aw 2 - 3(2 + e)A|uu| 

> (1 -e)X(u 2 + v 2 ), xeR 2 \D R£ . (3.46) 

Thus, using a comparison function argument and the property u 2 + v 2 = at 
infinity, we can obtain a constant C(e) > to make 



m 2 (x) + v 2 {x) < C(£)e-v / (i^l-l (3.47) 

valid. 

Let d denote any of the two partial derivatives, <9i and d 2 . Then (2.29) and 
(2.30) yields 

A(du) = \{2{du)e u - (dv)e v ), (3.48) 

A(dv) = 3X(dv)e v -2X(du)e u . (3.49) 

By computation and then using the Canchy inequality, we get 

A(\Vu\ 2 + \Vv\ 2 ) = 2[\V(d 1 u)\ 2 + \V(d 2 u)\ 2 + \V(d 1 v)\ 2 + \V(d 2 v)\ 2 } 

+4A|Vw|V + 6A|Vu|V - 2X8^0^ - 2Xd 2 ud 2 ve v 
-AXdmd^e* 1 - 4Xd 2 ud 2 ve u 
> 4A|Vu| V + 6A|Vu|V - A|Vw| 2 (l + 2e 2u ) 

-X\Vv\ 2 {2 + e 2v ), xeR 2 \D R . (3.50) 
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Therefore, as before, we conclude that for any e : < e < 1, there is a R £ > R, 
so that 

A(|V«| 2 + |Vt;| 2 ) > 2(l-e)(|V«| 2 + |Vv| 2 ), xeR 2 \D Rg . (3.51) 

Noting the property \Vu\ 2 + \Vv\ 2 = at infinity, applying the comparison prin- 
ciple, we arrive at 

|Vw| 2 + \Vv\ 2 < C{e)e~^ 2{1 ~ £)Mx \ \x\ > R. (3.52) 

Inserting this information into (3.4) and (3.5), we see that the associated functions 
u,w and the right-hand sides of (3.4) and (3.5) all lie in L 2 (R 2 ). Consequently, 
the pair of functions u and / yields a W^ 2,2 (lR 2 )-solution of (3.6) and (3.7), which 
must be the unique critical point of the functional I produced earlier. 
We may summarize our results as follows. 

Theorem 3.1 For any distribution of the points p±,P2, • • -p n G M?, the system of 
nonlinear elliptic equations (2.29) and (2.30) subject to the boundary condition 
(3.1) has a unique solution. Furthermore, the solution satisfy (3.47) and (3.52) 
decay estimates at infinity. 

4 Solution via variational approach on doubly 
periodic domain 



In this section, we consider solutions of (2.29) and (2.30) defined over a doubly 
periodic domain Q. In order to get rid of the singular source terms, we introduce 
a background function i> satisfying 

A n 

Av = -^+4nJ2S Ps (x)- (4-1) 

' ' s=l 

Using the new variable w so that v = i>o + w, we can modify (2.29) and (2.30) 
into 

Au = A(2e" - e va+w - 1), (4.2) 

4-7T77 

Au + Aw = 2X(e VQ+w - 1) + — . (4.3) 

Note that, since the singularity of v at p s is of the type In \x — p s \ 2 , the weight 
function e v ° is everywhere smooth. 
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To proceed further, we take u+w = f. Then the governing system of equations 
become 

Au = \(2e u -e V0+f ' u -l), (4.4) 
Af = 2X(e V0+f - u - 1) + (4.5) 

Integrating (4.5) and (4.4), we have 



/ 

Jn 



e V0+f - u dx = \n\ - = Ci > 0, (4.6) 

A 

/ e u dx = - [ e vo+f - u dx + -\n\ = + \n\) = C 2 >0, (4.7) 
Jn 2 J n 2 2 

Of course, the conditions (4.6) and (4.7) imply that the existence of an n- vortex 
solution requires that C\ > and C 2 > 0, which is simply 

\n\ - ^ = Cl > 0, (4.8) 
A 

since C\ > contains C 2 > 0. 

We can prove that (4.8) is in fact sufficient for existence as well. 

Theorem 4.1 The system of the non-Abelian vortex equations (4-4) an d (4-5) 
has a solution if and only if (4-8) holds or 

2vm < A|fi|. (4.9) 

Furthermore, if a solution exists, it must be unique, which can be constructed 
through solving a multiply constrained minimization problem. 

We use W 1,2 (Q) to denote the usual Sobolev space of doubly periodic functions 
over the cell domain Q. We will prove Theorem 4.1 in terms of three lemmas as 
follows: 

Lemma 4.1. Consider the constrained minimization problem 

min{/( W , f)\u, f e W^ 2 (Q), J k (u, f) = C k ,k = 1, 2}, (4.10) 

where 

I(u,f) = ^{l\Vu\ 2 + ±\Vf\ 2 -\u-\f+ 2 -^f}dx, (4.11) 

Ji(u,f) = [ e Vo e f - u dx = C 1 , (4.12) 
Jn 

MuJ) = [ e u dx = C 2 . (4.13) 
Jn 

Then a solution of (4-10) is a solution of the system of equations (4-4) an d (4-5)- 
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Proof. It is clear that the Frechet derivatives dJi, dJ 2 of the constraint function- 
al are linearly independent. 

Let (u, f) be a solution of (4.10). Then by standard elliptic regularity theory 
(u, f) must be smooth and there exist Lagrange multipliers Ai, A2 G K so that 

Au = -A- X 1 e Vo e f - u + A 2 e n , (4.14) 

4-7TT? 

Af = -2A + 2A 1 e^°e / " M + -p-. (4.15) 

Integrating the equation (4.15) and using Ji(w, /) = Ci, we obtain Ai = A which 
means that (u, f) verifies the equation (4.5). To recover the equation (4.4), we 
use J 2 (u, f) = C 2 . By virtue of Ai = A and integrating the equation (4.14), we 
have A2 = 2A. 

In particular, (it,/) is the solutions of the equations (4.4), (4.5). The lemma 
is proven. 

The admissible set of the variation problem (4.10) will be denoted by 

C = {u, f G W^MlMu, f) = C k , k — 1, 2}. (4.16) 

When (4.6) and (4.7) are satisfied, C U C 2 > 0. Thus C ^ 0. 

Lemma 4.2. If the condition (4-9) holds, then (4-10) has a solution. In other 
words, the system (2.29)-(2.30) has a solution if and only if (4-9) is fulfilled. 

Proof. By virtue of lemma 4.1, it is sufficient to show the existence of a minimizer 
of the constrained optimization problem (4.10). 

We first proved that under the condition (4.8) or (4.9), the objective functional 
I is bounded from below on C. For this purpose, we rewrite each 77 G W l,2 {VL) as 
follows 

r] = 77 + 7/, (4.17) 

where 77 denotes the integral mean of 77, 77 = J Q rjdx and f n rj'dx = 0. Hence, I 
may be put for (u, f) G C in the form 

I(u, f) = I {\\Vu'\ 2 + -\Vf'\ 2 }dx - A/|Q| + 2nnf - \u\Q\. (4.18) 
in 2 4 

Setting 

A(u,/) = -(X\n\-2nn)f-X\n\u, (4.19) 
we can derive from (4.6) and (4.7) the expressions 

u = lnC 2 -ln(y e M '), (4.20) 
/ = ]n(CiC 2 ) - ln (y e "') " ln (/ e "° +/ ' _U ')- ( 4 - 21 ) 
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Inserting (4.20) and (4.21) into (4.19), we have 

A(«,/) = (2A|H| -27m) In e u ') + (A|Q| - 2vm) ln^ e"° +/ '- u ') + C 3 , 



where C 3 = (2tto - A|fi|) InCi + (27m - 2A|ft|) lnC 2 . 
Using Jensen's inequality, we get 



In 



exp(v + /' - u) 



In 



Ye"'" > 



> In 
= In 
In 



|Q| exp 



vo + 



\Q \ exp 



\n\ 



u 



f - «')) 



in mi. 



Noting (4.9), we have 

A(«,/) > (2A|H| -27rn)ln|n| + (A|fi| - 27m) In 

+c 3 . 

Inserting (4.22) into (4.18), we arrive at the coercive lower estimate 



/(«, /) > j[ {^|Vu'| 2 + ^|V/'| 2 }dz - C 4 



(4.22) 



(4.23) 



where C4 > is an irrelevant constant. From (4.23), we know that the existence 
of solution of (4.10) follows. 

In fact, let {(uj, fj)} C C be a minimizing sequence of the variational problem 
(4.10) and set 



f ■ = T7TT / fjdx, Uj — T7TT / u jdx 



Q\ 



(4.24) 



Then, with u'j = Uj — Uj and f- = fj — f , we have u'j = and /'. = 0. In view of 



—j 



(4.23), we see that {(u'j,fj)} is bounded in W ' (Q). Without loss of generality, 
we may assume that {(u'j, /j)} converges weakly in W l,2 {Vt) to an element (uf, /') 
(say). The compact embedding 

W 1 ' 2 ^) ^ L P {Q), p>l, 



(4.25) 



then implies (u'j, /•) — > (u' , f) in L P (VL) (p > 1) as j — > 00. In particular, vf — 
and /' = 0. 

Recall the Trudinger-Moser inequality [1] 



(e F dx<C{e)ei^(\— + e]!\VF\ 2 dx\ FeW^ 2 (Q), F — 0. (4.26) 
In VL167T J J n J 
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where C(e) > is a constant. In view of (4.25) and (4.26), we see that the 
functionals defined by the right-hand side of (4.20) and (4.21) are continuous in 
u', f with respect to the weak topology of W 1,2 (Q). Therefore, Uj — > u, f — > f 
as j — > oo, as given in (4.20) and (4.21). In other words, (u, /) — (u + u', f + /') 
satisfies the constraints (4.6) and (4.7), and solves the constrained minimization 
problem (4.10). Thus Lemma 4.2 is proven. 

Now we state the uniqueness of the solution to the equations (4.4) and (4.5) 
as follows. 

Lemma 4.3. If system (4-4)-(4-5) has a solution, then the solution must be 
unique. 

Proof. Consider the following functional, 

1 1 27TT7 

J(u,f) = -\\Vu\\l + -\\Vf\\l + (-\u-\l+—m\ 

+ I {Xe V0+f - u + 2Xe u }dx. 
Jn 

It is straightforward to check by calculating the Hessian that J is strictly convex 
in W 1,2 (Q). Thus J has at most one critical point. However, any solution of 
(4.4)-(4.5) must be a critical point of J. This proves the lemma. 

5 Solution via fixed point method on doubly pe- 
riodic domain 

In this section, we shall solve the problem (2.29) and (2.30) by a fixed-point 
method via the Leray-Schauder theorem. This approach is of independent in- 
terest because the a priori estimates obtained in the process provide additional 
information on the governing equations. 

We rewrite (2.29) and (2.30) as 

Au = X(2e u -e v -l), (5.1) 

n 

Av = A(-2e n + 3e ,; -l)+4^^^ s (x). (5.2) 

s=l 

From Section 4, we know that (4.6) and (4.7) are a necessary condition for the 
solvability of (5.1) and (5.2) over a doubly periodic domain Q. 
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We now proceed to prove (4.6) and (4.7) are also sufficient for the existence 
of a solutions to the equations (5.1) and (5.2). Using a fixed point argument over 
the Sobolev space W 1,2 (Vt). Consider the proper subspace of W l ' 2 (Vt) defined by 

X = |( , u / ,f') G W 1,2 (Vt) J u'dx = 0, J v'dx = o}, (5.3) 

where u' — u — u, v' — v — v, and 

v = lnd-ln( 

u = In Ci — In | 



e v 'dx}, (5.4) 
e u 'dx^j, (5.5) 



after resolving the constraints (4.6) and (4.7). By the Poincare inequality [25], 
we may define the norm of X as follow 

\\(u',v')\\ x = ||V«'|| L2(n ) + llVu'H^^). (5.6) 

For each given (u',v') G X, consider the equations 

AU' = x(^--^--l), (5.7) 
AV > = x(=?^ + -¥^-l)+4*Y6 p (x). (5.8) 



s=l 



By (4.6) and (4.7), we see that the right-hand side of (5.7) and (5.8) have zero 
average value on Q. Therefore the equations (5.7) and (5.8) has a unique solution 
(U',V) G X. This correspondence, (u',v') i — > (U',V), gives us a well-defined 
operator T that maps X into itself, (U', V) = T(u', v'). 

Theorem 5.1 The system of equation (5.7) and (5.8) has a solution if and only 
if the conditions (4-6) and (4-7) are valid. 

We will prove Theorem 5.1 in terms of two lemmas as follows: 
Lemma 5.1. The operator T : X i — > X is completely continuous. 

Proof. Let (u' n ,v' n ) —¥ (u' ,v' Q ) weakly in X as n — > oo. Then (u' n ,v' n ) — > (u' ,v' ) 
strongly in D>(Sl) (p > 1). Set (U' n ,V$ = T(u' n ,v' n ) and ([«') = T(u' ,v' o ). 
Then 

^L e<dx L e<dx L eU ° dx L eV ° dxJ 

o) " x \j^ + j^x + j^r x j^cTx)- (5 - 10) 
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Multiplying (5.9) and (5.10) by U' n — U' and — Vq, and integrating by parts, 
respectively, we obtain 



/ 

Jn 

J 

Jn 



2C 2 e< 2C 2 e< 



Jn X In e<dx L eUodx 

Jn \j a e^dx J n e^dx 
3Cie< 3de v 'o 



In e<dx 



j^K-W*. (5-12) 

Note that the boundedness of {(u' n , v' n )} in X and the Trudinger-Moser inequality 
[1] imply that 

sup / e u '"dx < C < oo, (5.13) 
n Jn 

sup ! e<dx < C < oo. (5.14) 
n Jn 

Therefore, from (5.11), we obtain 

/ \V{U' n -U'^dx < ^{j^- ( e<\u' n -u' Q \\U' n -U' Q \dx 

^j^J/^-^K-u^dx} 

< x{^J^e<\u' n -u' Q \\K-U^\dx 

J^\v' n -v' \\U' n -U' \dx), (5.15) 



2d 



where u' n and v' n lie between u' n ,v' n and u' ,v' , respectively. In (5.15), we have 
used the inequalities 



J^e<dx>\n\e W Q^J^u' n dx) = 



and 



/ 

Jn 



f e<dx> |Q|expfi- f v' n dx) = 
n v l"l Jn ' 
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Applying the Canchy inequality and Holder inequality, and (5.13), we have 

/ e<\u' n -u' Q \\U' n -U' Q \dx < 1 / e 2 <\u' n -u' \ 2 dx+ £ - [ \K~Utfdx 
Jn Ae Jn a Jn 



+ ^1^(^-^)11^(0) 

< c e ||< - ^ ||! 4(Q) + ^f\\v(K - K)\\l Hn) . 



(5.16) 



Similarly, 



Jn 



- v'oWK - v \dx < c e \K - vXnw + ^f\\v(K - v )\\h iay (5-i7) 



Inserting (5.16) and (5.17) into (5.15), and letting e > be small enough, we 
have 

\MK ~ U^\\h ( n) < C(||< - ^oll| 4( ^) + \Wn - t/olli*(n)), (5-18) 

where C > is a constant. 
For (5.12), we have 

l|v(K - KOHi^ < c(||< - <||i 4(n) + IK - 4,111^)). (5.19) 

From (5.18) and (5.19), we arrive at 

IKK - 0*, K - vt)\\ x < c(||< - + IK - 4>lll^)), (5.20) 

where C > is a constant. This proves that (U' n ,V^) — > (Uq,Vq) strongly in X 
and the lemma follows. 

We now study the fixed point equation labeled by a parameter t, 

(u' t ,v' t )=tT(u' t ,v' t ), < t < 1. (5.21) 
Lemma 5.2. There is a constant C > independent oftE [0, 1] so that 

\\(u' t ,v' t )\\ x <C, 0<t<l. (5.22) 
Consequently, T has a fixed point in X . 
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Proof. When t > 0, it is straightforward to check that (u' t ,v' t ) satisfies the equa- 
tions 

Au' t = Xt(^^--^--l), (5.23) 

Av' t = \t(^^ + ^^-l)+4ntjrS Ps (x)- (5-24) 
Jn etdx Jn etdx jr[ 

In the doubly periodic domain Q, we let p, q G Q so that 

u' t {p) = max{u' t (x)\x G fl}, v' t (q) = max{v' t (x)\x G Q}. 

To facilitate our computation, we adopt the notation 

C e n * C e v>t 

J n e u tdx Jn etdx 

Then from (5.23), we have 

0>(Au[)(p) = Xt(2h' t (p)-g f t (p)-l). 

Therefore 

2/4(p) < fli(p) + 1 < + 1 = ^(«) + L 

Hence, for any i G fi, we have 

2^(x) < ^(qr) + 1, Vx G (5.26) 
From (5.24), using (5.26), we obtain 

g'M) < 1- (5-27) 
In view of (5.26) and (5.27), for any x G f2, we have 

^(x) < 1, h' t (x) < 1, x G ft. (5.28) 
Set t>£ = t^o + w' t . Then the equations (5.23) and (5.24) are modified into 

Au' t = At(-P^— - r \ — r- -1), (5.29) 

. , ^ ,-2C 2 e< 3Cie too+ ^ 4?m 
A«) = Xt(- r T-r+ r ' , , ~ 1 + -rprrt, 5.30 



where Av = + 47r^£ Pa (x). Multiplying (5.29) and (5.30) by u^iuj and 

s=l 

integrating by parts, respectively, and using (5.28), we have 
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||(V^,V^)lli 2( n) 



< 



xt 



de tv °- 



2C 2 e< 



f a e u 'tdx f n e tV0+w 'tdx 



Xt(- 



-2C 2 & 



+ 



3 Cie tv +w> t 



1) + 



dx 
Ann 



t\ ■ w'Adx 



J n e u 'tdx J n e tVt > +w 'tdx 
< I A(2 + 1 + l)\u' t \dx + I ((2 + 3 + l)A + -^K|cte 



< C\ \u' t \dx + C 2 / liuilda; 



< C £ + Cs\\(Vu f t: Vw f t )\\ 2 mn), 
Let e > be small enough, we have 

\\(u[,w[)\\ x = \\(Vu' t ,Vw' t )\\ LHn) <C, 



(5.31) 



(5.32) 



where C > is a constant. The existence of a fixed point is a consequence of 
Lemma 5.2, the apriori estimate (5.22) and the Leray-Schauder theory. The proof 
of the lemma thus follows. 



6 Further extensions 



In this section, we show that our method may be applied to establish the same 
existence and uniqueness theorem for multiple vortex solutions in the U(3) x U (3) 
model. 

We note that in the study [4] of the non-Abelian multiple vortex equations 
(2.17)-(2.20) the real- valued scalar field k and the complex- valued scalar field <fi 
are allowed to independently generate vortices with their respectively prescribed 
zero sets 

Z<p = {pi,P2,-,Pn}, Z K = {q 1 ,q 2 ,...,q m }. (6.1) 

In such a context, we can similarly develop an existence and uniqueness theory 
for the solutions of the equations by the same variational methods. To see this, 
we observe that, with the prescribed zero sets given in (6.1) for the fields k and 
and in terms of the variables u = Inn 2 and v = In |0| 2 , the governing system of 
nonlinear elliptic equations (2.29) and (2.30) is modified into 

m 

Au = \{2e u -e v - 1) + 4tt ^ ^(x), (6.2) 

m n 
t=l s=l 
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with the associated boundary condition 

u,v — >0 as \x\ — > oo. (6.4) 
Parallel to Theorem 3.1, we have 

Theorem 6.1 The system of nonlinear elliptic equations (6.2) and (6.3) subject 
to the boundary condition (6.4) has a unique solution for which the boundary 
condition (6.4) rnay be achieved exponentially fast. 

In order to prove Theorem 6.1, we introduce the background functions as before, 

m 

Mx) = -^hi(l + r|:r-g t |- 2 ), (6.5) 
t=i 

n 

v (x) = -^ln(l + r|x-p s |^ 2 ), r>0. (6.6) 



s=l 



Then 



Au = -hix) +4tt^^(x), (6.7) 

n 

Av = -h 2 (x)+4nJ2S P s(x), (6-8) 



s=l 



where 



hix) = 4V- — , (6.9) 

^ (r+ |x-%| 2 ) 2 ' 



t=i 



= *E (r+ |/- p , ffa - ( 61 »> 

We set -u = m + = + ^2, and f — wi + wi- Then (6.2) and (6.3) 
become 

Aiui = \(2e U0+Wl - e v « +f - Wl - 1) + h^x), (6.11) 

Af = 2\{e V0+f - Wl - l) +h 1 {x) +h 2 (x). (6.12) 

It can be checked that (6.11) and (6.12) are the Eular-Lagrange equations of 
the action functional 

I(w u f) = ^ jl-lVH 2 + TXI V /| 2 + 2(e M0+Wl - e u °) + (>+/-i - e"°) 
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It is clear that the functional / is C 1 over W 1,2 (M?) and strictly convex. We 
can use the methods in [16] and in the earlier study in the present paper to 
establish the coercive bounds 

DI(w 1 ,f)(w 1 ,f) > d ( || wi llwi^) + ||/|| w i,2 (M2) ^ -C 2 . (6.14) 

Therefore, it follows that the functional I has a unique critical point in W 1,2 (M. 2 ) 
which establishes the existence and uniqueness of a classical solution to the system 
of equations (6.2) and (6.3) subject to the boundary condition (6.4). 

We now turn our attention to the existence of multivortex solution over a 
doubly periodic domain Q. 

Take u and v over Q to satisfy 

Au = -^ + 47r^S 9t (x), &v = -^+47rJ25 Ps (x). (6.15) 

I I t=l I' s=l 

Then setting ^ = ^0 + ^15^ = ^0 + ^2, we see that the equations (6.2) and (6.3) 
over the doubly periodic domain Q become 

Airm 

A Wl = X(2e uo+wl - e V0+W2 - 1) + — , (6.16) 

Air 

A(w 1 + w 2 ) = 2X(e V0+W2 - 1) + — (m + n). (6.17) 

Theorem 6.2 For the vortex equations (6.16) and (6.17) defined over a doubly 
periodic domain fl, there is a solution if and only if the inequalities 

2n(m + n) < X\Q\, (6.18) 

7Y(3m + n) < A|Q|, (6.19) 

are satisfied. Moreover, if a solution exists, it must be unique. 

In the special case when the real scalar field k has no zero, that is, m = in 
(6.18) and (6.19), we recover Theorem 4.1. 

To proceed in the formalism of calculus of variations, we use the new variables 
g = w\, f — w± + w 2 . Then (6.16) and (6.17) take the form 

Ag = X(2e U0+9 - e Vo+f ~ 9 - 1) + (6.20) 

Att 

Af = 2X(e V0+f - 9 - 1) + — (m + n). (6.21) 

Integrating these two equations and simplifying the results, we arrive at the 
constraints 

r o-ir 
e V0+f ~ 9 dx = \n\ - — (m + n) = a ± , (6.22) 
Jn X 

j 4-7T777 

e U0+9 = -(a. + M-— -) = a 2 . (6.23) 

2 A 
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In order to show that the necessary condition ct\ > 0, a>2 > 0, which is exactly 
what stated in (6.18) and (6.19), is also sufficient for the existence of a solution, 
we recognize that Eqs. (6.20) and (6.21) are the Euler-Lagrange equations of the 
action functional 

I(f, 9 ) = ^{^|V/r + ^|V^| 2 + 2e« + e^ + ^ 

{Anm \ /27r(m + n) \ A , . 

+ (w~ 1 ) 9+ ( _ W~ 1 ) / ' <te - (6 ' 24) 

Now decompose /, g into f — f' + f,g — g' + g with f,g<E~M, and f Q f'dx = 
Q,J n g'dx = 0. Thus, applying (6.22) and (6.23), we may rewrite (6.24) in the 
form 

KM - /^iv/f + ^ivsf}^ 

= ailn(/" e V0+f '- 9 'dx) + 2a 2 ln( f e Uo+9 ' dx) 
Jn Jq 

+«i(1 -In ai) + 2a 2 (l -In a 2 ). (6.25) 

It is seen immediately that the right-hand side of (6.25) has a uniform lower 
bound in view of the Jensen inequality again. So the existence of a critical 
point of (6.24) subject to the constraints (6.22) and (6.23) follows as before. The 
uniqueness of a critical point of (6.24) results from the convexity of the functional. 
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